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We discuss how variations in the scanning tunneling microscope (STM) tip, whether unintentional 
or intentional, can lead to changes in topographic images and dl/dV spectra. We consider the 
possibility of utilizing functionalized tips in order to improve the sensitivity of STM experiments to 
local irregularities at the surface or hidden below the surface layers. The change in the tip symmetry 
can radically alter the contrast of the topographic image due to changes in tip-surface overlap. The 
dl/dV curves change their shape according to which sample bands the tip orbital tends to overlap. 
In addition, relative phases between competing tunneling channels can be inverted by changing the 
tip symmetry, which could help reveal the origin of a local irregularity in tunneling spectrum. 



PACS numbers: 68.37.Ef 71.20.-b 74.50. +r 74.72.-h 

I. INTRODUCTION 

Scanning tunneling microscopy and spectroscopy 
(STM/STS) are extensively used to probe the quasipar- 
ticle spectra of high temperature superconductors, espe- 
cially of the cuprate materials such as Bi2Sr2CaCu20s+^ 
(Bi2212)PH^. These studies have contributed immensely 
in obtaining insight into superconducting and pseudogap 
phases of these materials owing to the very high spatial 
and energy resolution of the spectroscopy. The full ex- 
perimental potential of STM/STS is, however, far from 
exploited, since the chemical and spatial resolution as 
well as the directional selectivity could, in principle, be 
improved by overcoming various challenges related to in- 
trinsic properties of STM. 

The first challenge is related to the fact that the sig- 
nal from superconducting layers is filtered by insulating 
oxide layers. In addition, the signal does not arrive at 
the microscope tip directly from the site below the tip 
but rather through symmetry driven tunneling channels. 
Hence, experimentally observed spectral features are not 
only distorted by filtering, but ascertaining their spatial 
origin becomes harder the deeper within the sample is 
their source. Another factor to be aware of is that the 
electronic and geometrical properties of the microscope 
tip may also affect STM images. In standard interpreta- 
tion, the tip is assumed symmetric to rotations around 
the vertical axis ('s-wave'). The relative blunt ness of 
the tip suggests this assumption and, furthermore, the 
metallic nature of common tip materials does not favor 
oriented "dangling" orbit als, and specific effects of tip 
geometry and electronic structure are usually ignored. 
However, this interpretation has been challenged by a 



variety of sophisticated theoretical studied^', which, on 
the one hand show that the closer the tip is to the sur- 
face, the more important the symmetry of tip orbitals is 
to the STM image and, on the other hand, they indicate 
that s-wave tips might not be able to produce images 
with atomic resolution. Notably, the relative importance 
of tip orbitals with different symmetries in the STM sig- 
nal can be controlled by varying the tip-surface distance, 
or alternatively, by functionalizing the tip by attaching 
particular atoms/molecules at the apex of the tip. 

In various experimental studies^^, the tip proper- 
ties have been deliberately modified for better spatial or 
chemical selectivity by picking up a molecule, e.g. CO, 
O2 or N2, from the surface to the apex of the tip af- 
ter which the plain tip results have been compared to 
the results with a functionalized tip. This has proven 
to improve chemical selectivity in adsorbate on a solid 
surface system since each surface species responds in an 
individual way to the change of tip symmetry due to 
functionalization. Even in the case of relatively compli- 
cated molecules, such as pentacene attached to the tip, 
the contrast of differential conductance maps revealing 
the geometric shape of molecular orbitals is surprisingly 
sharpened^ by the enhanced directionality of the probe- 
sample overlap due to the dominance of the a or 7r type 
localized orbitals closest to the sample surface. The idea 
of utilizing a functionalized tip is further encouraged by 
the fact that the resolution of atomic force microscopy 
(AFM) studies is enhanced by using CO molecule at- 
tached to the metallic tipP^, although improving the res- 
olution and selectivity in the case of AFM has a different 
physical basis from STM. 

The preceding work on molecular adsorbate systems 
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provides motivation for considering functionalized tips 
for investigating cuprate materials, where the interesting 
physics is found at layers not immediately on the surface. 
We can foresee many uses of functionalized tips in prob- 
ing these materials. An obvious application would be a 
more accurate probe for the surface properties. In the 
case of Bi2212, to date one has mostly observed the Bi 
atoms of the surface layer in topographic STM images. 
Information from the oxygen atoms of the surface layer 
has been difficult to obtain, and hence STM provides 
limited ability to observe the real symmetry of the 
surface. A more involved application is related to the 
dl/dV spectra. As shown in our previous studied, 
the spectra are mappings of the electronic structure of 
the cuprate layer, filtered by the BiO and SrO layers, 
which sit on top of the cuprate layers. Although the 
signal can be theoretically resolved into partial signals of 
different origin, the STS experiment measures a signal 
combined from different sources and, furthermore, the 
information from different directions in momentum space 
is integrated into a single outgoing signal. A possible 
means of decomposing the STS spectrum experimentally 
is to measure topographic spectral maps. In practice, 
real space and Fourier transformed maps differentiating 
electron and hole states have been produced in order 
to recognize direction dependent underlying electronic 
order 20 . A third possible application would be inves- 
tigation of symmetry properties related to nematicity 
transitions in the subsurface electronic structure at low 
values of doping, where the oxygen atoms of Cu02 layers 
in different directions appear non-equivalent below the 
pseudogap transition temperature T* as shown in Ref. 



In fact, functionalized tips may well have been used 
inadvertently in studies of Bi2212. It has been reported 
that tip resolution can improve dramatically after use, 
and it has been suggested that this may be due to the 
attachment of a Bi atom to the tunneling tipP^. Here 
we will explore the effects of tips, which are dominated 
by a single atomic orbital of varying symmetry. We find 
that such tips can not only couple with different surface 
orbitals, but can also affect the interference between tun- 
neling channels from the cuprate layer to the microscope 
tip. This suggests possibilities for investigating irregu- 
larities within the cuprate layers, which are not directly 
accessible to the STM tip, and thus reveal features invis- 
ible to STM/STS experiments with a symmetric tip. In 
this way, our study advances the understanding of matrix 
element effects in STMpM] 

This article is organized as follows. Section II briefly 
describes the model utilized in this study. Section III 
presents calculations of the effect of tip symmetry on to- 
pographical maps, STM corrugations and dl/dV spectra 
at various high symmetry sites. In Section IV we show 
how a hidden irregularity within the cuprate layer can 
be studied by comparing results obtained with different 
tip symmetries, among other results. Finally, conclusions 



are given in Section V. 



II. DESCRIPTION OF THE MODEL 

Although the model used is essentially the same as in 
Ref. 19, we provide an overview here for the sake of 
completeness. The simulation supercell consists of 8 unit 
cells that have been rotated 45 degrees as shown in Fig. 
[I] The unit cell geometry is adapted from Ref. [28] In 
total there are 120 atoms and 2 x 464 orbitals accounting 
for both spins. The topmost 7 layers below the cleavage 
plane are included, which are in descending order: BiO, 
SrO, CuOi, Ca, CuO^, SrO and BiO. Depending on 
the atom type different sets of orbitals are included. For 
Bi, O and Ca the orbital set is (s,p x ,p y ,p z ) and for Cu 
the orbital set is (4s, d^ z i_ r i, d xy , d xz , d yz , d x 2_ y 2). Sr 
has a single s orbital. The Green's function is computed 
in a supercell Brillouin zone at 256 equally distributed 
k-points corresponding to 8 x 256 = 2048 k-points in the 
unit cell. 

As shown in Fig. [TJ we define the position of Cu (or 
surface Bi) atoms as top site (T), the position of bonding 
O atoms of the cuprate layer as bridge site (B), and the 
position of O(Bi) as the hollow site (H). In this choice 
of geometry, the p x and p y orbitals point along the bond 
between the neighbouring Cu atoms, i.e., in the direction 
from the top site to the bridge site (T — » B), whereas the 
linear combinations ~^(Px ^Py) point from the top site 
towards the hollow site (T — » H). Concerning the oxygen 
atoms related to the cuprate layer, the apical oxygen is 
at the top site, while the bonding oxygens between Cu 
atoms are at bridge sites. 

The tip is modeled as a single orbital in order to high- 
light the effect of each type of tip symmetry on the 
STM/STS spectrum. We consider four different tip sym- 
metries which can be divided into two groups. These 
consist of s- and p^-orbitals, which are symmetric to ro- 
tation around z-axis, accompanied with two choices of 
horizontal p-orbitals. We use a shorthand notation p^ for 
a p-orbital along the bond between two neighbouring Cu 
atoms, and p r for an orbital rotated by 45°, i.e., aligned 
in the direction connecting next nearest neighbouring Cu 
atoms. These two horizontal tip symmetries can also be 
connected to the nodal (p r ) and antinodal (pt,) directions. 

The Hamiltonian is the same as that use d pre viously in 
our STM studies of optimally doped B22l3±^MH, w here 
superconductivity is modeled with a pairing matrix A a p 
between electrons and holes with opposite spins: 



(1) 



with real-space creation (annihilation operators) 
(c ao -). Indices are combination of site and spin indices. 
The d-wave symmetry of the superconducting gap is 
modeled by taking A to be non-zero only between the 
neighbouring d x i_ y i orbitals, and with opposite signs for 
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FIG. 1. (color online) (a) Geometry of the simulated system. 
Left hand panel gives the arrangement of the CuCb layer, 
while the right panel depicts the BiO and SrO layers on top of 
the CuC>2 layer; (b)-(e) Topographic STM images for different 
tip symmetries, (b) s-orbital, (c) p^-orbital, (d) p r -orbital, 
(e) pb-orbital. Indicated also are Bi and O positions and the 
top (T), bridge (B) and hollow (H) sites. In (d) and (e) the 
tip symmetry is further indicated with a schematic of the 
corresponding p-orbital. 



CherPand a similar approach of Sacks and NoguersPwho 
utilize derivatives of the real space tip wave function. 
There is an alternative way of writing Eq. (J2|: 
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J2 Pu\E F )M t , cPcd [E F + eV)M\, v (3) 



where 



M tc = V ts G°+V fc . 



(4) 



The difference between the two formulae is that Eq. ([2| 
emphasizes the coupling between the tip and the surface 
orbitals expressed by the hopping integrals V ts: whereas 
Eq. ([3| makes more transparent how tunneling takes 
place through the filtering layers between the tip and the 
superconducting layers. Of course the symmetry of the 
tip itself is described by V tsi but in order to recognize 
the relevant surface orbitals, it is essential to consider 
the filtering matrix elements, M tci as well. 

The spectral function p cc > depends on both regular and 
anomalous Green functions as shown below 
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To conclude, the formulae given above provide not only 
the total dl/dV spectrum, but also an explicit form of 
the filtering effect on the cuprate layer, and decompose 
the density matrix of the cuprate orbitals into regular 
and anomalous spectral terms. 



III. RESULTS 



A. Topographic maps 



x- and ^/-directions, which leads to a d-wave symmet- 
ric superconducting gap 32 . The tight-binding parameters 
V a p are taken to be doping-independent in the spirit of 
the rigid band model P^™^ 

The tunneling current is computed within the Todorov- 
Pendry approximation^* 4 ^ as 

a = % = £ P*(E F )V t , sPs AE F +eV)V s \ t . (2) 

tt'ss' 

The hopping integrals V a p depend on the symmetry of 
tip and surface orbitals as well as the distances involved 
through the Slater-Koster coemcient d 42 | 43 i The tip to 
surface distance is usually greater than the distance be- 
tween the atoms in the bulk material. To allow for this 
the standard Slater-Koster method is enhanced by an ex- 
ponential decay at large distances^. The Slater-Koster 
method allows one to take the directionality of both the 
tip and surface orbitals into account. It can be considered 
as a special case of the more general "derivative rule" of 



The topographic maps of Fig. [l[b-e) are calculated 
using Eq. Q in constant current mode, i.e., there is a 
feedback loop varying the tip-surface distance in order to 
retain the tunneling current constant. For bias voltage 
we choose —0.3V. We also follow the common practice of 
accompanying topographic maps with corrugation curves 
in Fig. [2j since the latter reveal the absolute scale of 
the contrast otherwise obtained in a rather arbitrary way 
from topographic maps. 

Experimental topographic maps typically resemble the 
calculated maps for s-wave tip with a high intensity at 
the position of the surface Bi atoms. Within the Tersoff- 
Hamann type analysis 40 , a high intensity signal is at- 
tributed to a high local density of electrons. As Fig. 
|3ja) indicates, the surface Bi atoms have a large density 
of states over a wide range of energies around the Fermi 
energy. This was demonstrated already in our earlier 
workP2 where we further related this to the coupling be- 
tween the p z orbitals of Bi and the d z i orbital of the Cu 
atom below through the intermediate p z orbital of the 
apical oxygen. In contrast, orbitals of O(Bi) contribute 
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FIG. 2. (color online) Corrugation curves for various tip or- 
bitals along the route (A — >• B — >• C —> A) shown in the 
insert. The vertical scale represents the height change the tip 
experiences along the route. Curves for p r and pb are similar 
in overall shape yet they differ strongly in relative amplitude. 
This shows that they connect at least partially to different 
surface orbit als. Note that each curve was shifted downwards 
by its own average tip surface distance thus making all curves 
zero centered. 



Fig. [2] Interestingly, the horizontal tip orbitals not only 
invert the contrast of the topographic image, but they 
are also more sensitive to the horizontal position than 
tips with rotational symmetry around the z-axis. This 
is seen in the more strongly varying corrugation as the 
tip is scanned between the top, bridge and hollow sites. 
These differences, originating from the symmetry depen- 
dent nature of the tip-surface overlap, might be useful in 
pursuing a better lateral resolution in STM imaging. 

Unfortunately, our theoretical predictions cannot be 
compared directly with measurements at this time since 
we are not aware of any experimental STM literature on 
Bi2212 where a deliberately functionalized tip has been 
deployed. However, some insight can be obtained by con- 
sidering a study such as that of Ref. [T7J which discusses 
the effects of attaching O2 functionalizing molecules to 
a tungsten tip for investigating a metal phtalocyanine 
monolayer on an Au(lll) surface. Ref. [TT] finds selec- 
tivity with respect to molecular orbitals between normal 
and 02-functionalized tips, and observes the STM pat- 
tern of the adsorbate monolayer to rotate as different tips 
are used to probe the surface. 



to the BiO bands, which contain a gap around the Fermi- 
energy. In Fig. J3^a) these states are broadened since we 
have invoked a Fermi liquid type self-energy in the non- 
cuprate layers^. 

As seen in Figs. [l](b) and[l](c), tips with s- as well as 
p^-symmetry expose the Bi atoms as bright spots. This 
is hardly surprising, especially in the case of the rather 
strongly oriented p z tip orbital. Our computations also 
suggest that p z symmetry slightly sharpens the contrast 
of the topographic map, but as shown in the corrugation 
curves of Fig. [2j the effect may not be experimentally 
significant. Tips with in-plane p-symmetry show striking 
changes in topographical maps. In the case where the 
tip orbital lies along the top-hollow (nodal) direction, p r 
(Fig. [ljc)), the contrast is inverted so that O(Bi) atoms 
appear as bright spots. However, as discussed further in 
Section IV below, visibility of oxygen atoms is spurious 
since the tunneling signal is still dominated by overlap 
between the horizontal tip orbital and the p z orbitals of 
Bi atoms. This is also the obvious reason why the bright 
spots in the topographical image lose their sharpness. 
The change in contrast becomes even more dramatic in 
the case of the tip orbital lying along the top-bridge (anti- 
nodal) direction, p^. As seen in Fig. |l](d), the fourfold 
symmetry is now broken, and bright lines appear per- 
pendicular to the orbital direction. In addition, intensity 
maxima are seen at bridge sites when the bonding oxy- 
gen orbital points to the same direction as the tip orbital. 
This indicates that in this case also the nature of the to- 
pographic image is driven the overlap between the tip 
orbital and the p z orbital of Bi atoms. 

The intensity patterns in topographic images can be 
quantified as corrugation curves, such as those shown in 



B. dl/dV spectra 

The shape and intensity of a dl/dV spectrum depend 
on the tip symmetry as well as its position. In Figs. [3jb- 
d), the dl/dV curves are calculated at 5 A vertical dis- 
tance between the tip and the surface atoms. At this tip- 
surface distance, the gap features are clearly observed, 
and calculations indicate that the signal predominantly 
comes through the p z orbitals of the Bi atoms. There are 
clear changes in the VHS features at ~-0.15 to -0.4 eV 
with tip site and orbital symmetry: they seem to vanish 
for the s-wave tip at a hollow site, and for both the hor- 
izontal p-wave tip orbitals at a top site. The apparent 
reason for this is the increasing spectral weight of the 
BiO bands with increasing binding energy. These config- 
urations seem to favour a relatively large overlap between 
the tip orbital and the p-orbitals of O(Bi) atoms. For an 
s- wave tip on the top site, the dl/dV curve resembles the 
Bi p z LDOS, as expected. However, this is also true for 
the horizontal p-wave tips at the bridge or hollow site. 
These results support the conventional assumption that 
the experimental tips are predominantly s-wave. 

We note also a trend in the overall intensities of the 
spectra. As expected, horizontal p-symmetry of the tip 
orbital tends to decrease spectral intensity due to de- 
creasing overlap. For s-wave symmetry, the intensity is 
highest for the top site and the lowest for the hollow site. 
For the horizontal p-symmetry, the top site gives the low- 
est intensity due to antisymmetry between p z of Bi and 
the tip-orbital. Moving the tip either to bridge or hollow 
site allows an increased overlap between the tip orbital 
and p z of Bi. 
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FIG. 3. (color online) (a) LDOS ofp* of Bi, p x/y of O(Bi) and 
p z of O(Bi). The first orbital hybridizes with the C11O2 bands, 
the second is dominant in BiO bands (Bi pocket states), and 
the third is rather localized; (b), (c) and (d) show the dl/dV 
curves with different tip orbitals at the high symmetry sites 
B, H and T, respectively. 



C. Impurity mapping and subsurface irregularities 

dl/dV spectra are also used to detect and analyze hid- 
den subsurface irregularities. These can include sub- 
stitutional impurities, such as Zn or Ni replacing a Cu 
atonP, or dopant atoms like oxygen, the fingerprints 
of which have been discussed both experimentally 2 ^ and 
theoreticalljEl However, perhaps the most fascinating 
potential application may be the C4 —> C2 nematic in- 
tracell symmetry breaking of the subsurface electronic 
structure observed at very low values of hole doping in 
the pseudogap phased This phase transition in the elec- 
tronic order is not observed in the topographic images, 
but rather in dl/dV maps and in the so called Z-maps 
derived from the differential conductance data. The rel- 
evance to the present theoretical considerations arises 
from the fact that the subsurface bonding oxygen atoms 
in Cu02 appear non-equivalent when this low-symmetry 
phase is mapped. We propose that breaking the symme- 
try of the tip by functionalization would give additional 
information about the nematic phase. Note that the ob- 
served patterns are difficult to interpret due to the fil- 
tering effect of the BiO and SrO layers. It is interesting 
therefore to explore whether different tip symmetries can 
provide insight by probing different tunneling channels. 
In earlier work 19 we showed that the bonding oxygens of 
the Cu02 layers are invisible to a spherically symmetric 
tip placed right above that atom. Hence we choose a local 
perturbation at this bonding oxygen. The simplest kind 
of perturbation is a change of the onsite energy, which 
could be induced by a be nearby dopant atom. 

Figure (4^d) shows a tunneling tip positioned above a 



bonding oxygen in the cuprate layer, equivalent to the 
bridge site of Fig. 1. To simulate the effect of a nearby 
impurity, we modify the onsite energy of the p-orbital 
along the bond by -2.0 eV. 4 ^In Fig. Etc) we compare its 
LDOS (blue line) to that of the regular site (red line). 
The dl/dV spectra from the bridge site with and with- 
out onsite modification are plotted for s (Fig. [4] a) and 
Pb Qb)) type tips to demonstrate how the impurity ef- 
fect modifies the spectrum in different tunneling chan- 
nels. Although the difference between the two LDOS 
results is striking (Fig. [4^c)), the spectra obtained by a 
rotationally symmetric tip do not reflect this difference. 
This can be explained by the results of Ref. HH where 
we showed that two tunneling channels in antiphase are 
open between the oxygen and the tip (Fig. |4jd)), causing 
destructive interference between the corresponding sig- 
nals. [The observed signal comes from different channels 
involving the Cu atoms.] 

In Fig. (4^d) we show the dl/dV curve measured by a 
tip with pb orbital, and a clear difference is seen between 
the spectra of the regular and perturbed oxygens. The 
shapes of LDOS (b) and the spectra (d) are not identical, 
but the relative change of the spectrum follows that of 
LDOS, especially over the range between —0.5 — O.OeV. 
The change in the spectrum of Fig. |4jd) is caused by a 
reversal of the relative phase between the two tunneling 
channels in Fig. |4ja) due to p-type symmetry of the 
tip, a point to which we will return below. This result 
suggests that experiments probing the breakdown of C4 
symmetry, such as those in Ref. 21, could be further 
clarified by varying the tip symmetry. 



IV. ANALYSIS 

A. Tip symmetry and contrast changes in the STM 

map 

A proper understanding of the observed topographical 
images and spectra requires an analysis of the filtering 
matrix elements M tc which describe the coupling between 
the tip and the surface orbitals, propagation across the 
oxide layers, and coupling between the superconducting 
cuprate layer and the covering oxide layers. Starting from 
the latter, the cuprate layer is coupled to the oxide lay- 
ers mainly through the overlap between the d z i orbitals 
of the Cu atoms and the p z orbitals of the apical oxy- 
gen above. The propagation through the oxide layers 
is mediated by p z orbitals of the apical oxygen and the 
Bi atom above on the surface. When the tip symme- 
try is varied, the relevant factor is that p z of the surface 
Bi is the dominant surface orbital mediating the signal 
from the cuprate layer to the microscope tip. Hence, 
the analysis of the observed results can be built on the 
overlap between the tip orbital and the p z orbitals of 
the Bi atoms in the vicinity of the tip. While other sur- 
face orbitals are present, they are not strongly coupled 
to the cuprate bands, and hence they have a relatively 
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low spectral weight in the vicinity of the Fermi energy. 
Thus, the horizontal p-orbitals of Bi and O(Bi) atoms 
are the dominant orbitals of the BiO bands related to the 
"bismuth pocket" states above Fermi-level and the corre- 
sponding valence bands below the "spaghetti" of bands 
(Fig.ga)^l. 

Since only the Bi p z surface orbital is important, the 
analysis of the various tip orbitals is rather straightfor- 
ward. The s— and ^-orbitals couple best on the top site, 
as seen in Fig. [3^b). For the horizontal tip orbitals, the 
overlap between the tip and p z of a surface Bi is lowest 
when the tip is directly above the Bi atom. In the case of 
p r both the bridge and hollow sites seem to give a large 
overlap between tip orbital and p z of Bi, but the hollow 
sites appear brighter since the overlap between the hor- 
izontal p-orbitals contributing to BiO bands is larger at 
this site. For pb the fourfold symmetry is broken. Since 
the orbital is oriented along the bonds between Cu atoms, 
there is a zero overlap line when following the bonding di- 
rection perpendicular to the orbital orientation. On the 
other hand, pb overlaps strongly with the two Bi atoms 
at the bridge site in the antinodal direction of pb. Mov- 
ing towards the vicinal hollow sites preserves the strong 
overlap. At the hollow site, a strong overlap with four 
Bi atoms can be found. Since the distance is V^-fold as 
compared to the bridge site, overall the overlap remains 
roughly the same when moving from the bridge site to 
the hollow site, resulting in bright and dark stripes in 
Fig. [3e). 



B. Hidden irregularities and the relative phase 
between two tunneling channels. 
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FIG. 4. (color online) (a) dl/dV spectra of the regular and 
modified bonding oxygen using s-symmetric tip above the 
oxygen; (b) The same spectra using p^-symmetric tip above 
the oxygen; (c) LDOS of the p-orbital of the bonding oxygen 
in the cuprate layer. Regular case (red) is compared with an 
orbital with a modified onsite energy (blue); (d) Two compet- 
ing channels from the bonding oxygen (o) through d z orbitals 
(a or 6), apical oxygens (/ and g) and the surface bismuths 
(1 and 2) to the tip (t). 



and G b, and hence the two first spectral terms are equal 
to each other, but the last one has the opposite sign. 

Finally, using equations ([3]), Q and ([5| we obtain the 
total current from orbital o to the tip: 



We consider the possibility of analyzing subsurface ir- 
regularities by functionalizing the tip. This may be pos- 
sible if the interference between two different tunneling 
channels can be controlled. As an example, let us con- 
sider the situation where the tip is at the bridge site (see 
Fig. [4|(d)), i.e., above an oxygen atom in the Cu02 layer 
at an equal distance from the two neighbouring Bi atoms 
(1 and 2). This opens two channels from the Cu02 layer 
(d Z 2 orbitals of Cu atoms a and b) with filtering matrix 
elements M la = V t iG lf V fa and M 2b = V t2 G 2 ) g V g b which 
are equal apart from a phase factor. It is obvious that 
Mi a and M 2 b are equal in the case of s- or ^-symmetric 
tip, since in that case V t \ = V t 2, but the opposite sign is 
obtained with horizontal p-orbitals. 

To consider interference effects, we also need to com- 
pare elements of the spectral function 

q _ r<+ n 

u aoa — ^ao^oo^oa 
Obob — ^bo^oo^ob 
Saob ~ ^ao^oo^ob^ 

where o denotes the p-orbital of the oxygen along the 
bond between the Cu atoms. It is evident that the p- 
symmetry leads to an antiphase between elements G ao 



h =Ptt(\M la \ 2 S aoa + \M 2b \ 2 S b ob+ 

M la S aoh Ml 2 + MibSboaMl) (6) 
=p tt (2\M la \ 2 T2\M la \ 2 )S aoa . 

It is clear that '-' sign on the right hand side of Eq. 
(6) applies to the case with Vti = V t 2, i.e., the s- or 

symmetry. Hence the oxygen below the tip would be 
invisible in STM. On the other hand, constructive in- 
terference is obtained for horizontally p- symmetric tip. 
This should emphasize any variations in the horizontal p- 
orbital of the bonding oxygen. The situation is of course 
reversed if the bonding oxygen is replaced by an s-wave 
symmetric perturbation, i.e. there is a constructive inter- 
ference for s- and p 2 -symmetric tips. This explains our 
comparison of a spectrum at the bridge site for unper- 
turbed system with a vacancy at the site of the bonding 
oxygen. 

With regard to the C4 — »■ C2 transition, the main prac- 
tical conclusion of our analysis is that obtaining informa- 
tion related to the non-equivalent bonding oxygens can 
be enhanced by manipulating the tip symmetry by func- 
tionalizing molecules. This suggests that an s-wave sym- 
metric tip gathers the signal from a neighbourhood of 
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the bonding oxygen below, while a p-wave tip should be 
more selective to the on-site oxygen orbitals. 

V. CONCLUSIONS 

We have investigated the effect of different tip sym- 
metries on topographic STM images and dl/dV spectra. 
The change in tip symmetry can lead to contrast inver- 
sion and even to the breaking up of the C4 symmetry 
due to strong differences in the nature of the tip-surface 
overlap. Furthermore, the corrugation curves show that 
the horizontal p-orbitals at the apex of the tip would 
be more sensitive to position of the surface atoms. The 
dl/dV curves reveal the leading sample states (C11O2- 



band, BiO-band or some more localized states) to which 
different symmetries tend to couple. More strikingly, 
changing the symmetry of the tip can be used to change 
the relative phase between two open tunneling channels. 
This could be applied to detect subsurface irregularities 
which would otherwise be invisible to STM. 
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